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ABSTRACT 

We report on an X-ray flare detected on the active binary system II Pegasi 
To . with the Swift telescope. The event triggered the Burst Alert Telescope in the 

' hard X-ray band on December 16, 2005 at 11:21:52 UT with a 10-200 keV lu- 



minosity of 2.2x10^^ erg s~^ — a superfiare, by comparison with energies of 
typical stellar flares on active binary systems. The trigger spectrum indicates 
^ a hot thermal plasma with T~180 xlO^K. X-ray spectral analysis from 0.8-200 

keV with the X-Ray Telescope and BAT in the next two orbits reveals evidence 
for a thermal component (T>80 xlO^K) and Fe K 6.4 keV emission. A tail of 
emission out to 200 keV can be fit with either an extremely high temperature 
thermal plasma (T~3xlO^K) or power-law emission. Based on analogies with 
solar flares, we attribute the excess continuum emission to nonthermal thick- 
target bremsstrahlung emission from a population of accelerated electrons. We 
estimate the radiated energy from 0.01-200 keV to be ~6xl0^^ erg, the total ra- 
diated energy over all wavelengths ~10^* erg, the energy in nonthermal electrons 
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above 20 keV ~3xl0^° erg, and conducted energy <5xl0^^ erg. The nonthermal 
interpretation gives a reasonable value for the total energy in electrons > 20 keV 
when compared to the upper and lower bounds on the thermal energy content of 
the flare. This marks the first occasion in which evidence exists for nonthermal 
hard X-ray emission from a stellar flare. We investigate the emission mechanism 
responsible for producing the 6.4 keV feature, and find that coUisional ionization 
from nonthermal electrons appears to be more plausible than the photoionization 
mechanism usually invoked on the Sun and pre-main sequence stars. 

Subject headings: stars: activity — stars: coronae — stars: flare — stars: late- 
type — stars: flare — stars: individual (II Peg) — X-rays: stars 

1. Introduction 

The general model for solar (and stellar) flares involves a release of free energy via mag- 
netic reconnection, which accelerates particles and causes subsequent plasma motions and 
heating (Dennis & Schwartz 1989). Particle acceleration during solar flares can be diagnosed 
directly through nonthermal hard X-ray bremsstrahlung emission and radio gyrosynchrotron 
emission. Until now, evidence for particle acceleration during stellar flares has been either 
indirect (through proxies such as optical and UV flare emissions; Hawley et al. 1995, 2003) or 
achieved through the detection of nonthermal radio flare emissions (Klein & Chiuderi-Drago 
1987), under the assumption that the observed radiation is gyrosynchrotron emission from 
accelerated particles. The interpretation of radio observations, while useful for constraining 
the characteristic magnetic field in the source, is complicated by optical depth effects and 
spatial inhomogeneities in the emitting source, which are usually not uniquely determined. 
Previous detections of hard X-ray emission from stellar flares could not demonstrate unequiv- 
ocally the presence of nonthermal emission; instead, a supcrhot thermal component could 
equally well explain the observed 10-50 kcV emission (Favata & Schmitt 1999; Pallavicini 
2001; Franciosini et al. 2001) without the need for any additional power-law components. 
The detection and characterization of nonthermal hard X-ray emission during stellar flares 
would be important, as it would allow a direct investigation of the energy spectrum of acceler- 
ated electrons, and a comparison of the total thermal/nonthermal energy budget. Prospects 
for detecting such emissions have normally been taken to be remote, as during solar flares 
the nonthermal hard X-ray emission (20 keV- 1 MeV) is >10~^ times less intense than soft 
X-ray (1-10 kcV) emission (see, e.g., the composite solar flare spectrum in Fig. 94 of As- 
chwanden 2002); this is a consequence of coUisional energy losses exceeding bremsstrahlung 
losses by a factor of 10^. This necessitates large flare sizes for hard X-ray detections, which 
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are intrinsically rare, which coupled with a dearth of pointed observations exacerbates the 
situation. 

Iron K fluorescence emission has been seen during solar flares (Parmar et al. 1984; Zarro 
et al. 1992) and also seen in impulsive flares from some pre-main sequence stars (Tsujimoto 
et al. 2005). The origin on the Sun is attributed to photoexcitation of K shell electrons in 
photospheric iron by X-ray bremsstrahlung radiation. The emission mechanism for pre-main 

sequence stars is the same as for the Sun, but the fluorescence is presumed to arise instead 
from X-ray irradiation of an accretion disk rather than the stellar photosphere. There have 
not been any previous reports of iron Ka fluorescent emission from non-accreting stellar flare 
sources. Such an identification would provide another link between solar fiare physics and 
the physics of flares on other stars. 

II Peg (HD 224085) is a well-known active binary system which displays evidence for 
starspots covering a large fraction (>40%; Marino et al. 1999) of the primary's photosphere, 
and has signatures associated with vigorous magnetic activity across the entire electromag- 
netic spectrum. At a distance of 42 pc (Perryman et al. 1997), this single-lined spectroscopic 
binary is composed of a 0.8 Mq K2IV primary tidally locked in a 6.7 day orbit with a ~O.4M0 
secondary (Berdyugina et al. 1998). Space velocities indicate that this binary is a member 
of the old disk stellar population (Eggcn 1978). Both photospheric and coronal metal abun- 
dances have been determined, leading to an apparent depletion of coronal iron by about 
a factor of 4 relative to the photospheric value (Huenemoerder et al. 2001, and references 
therein) ^. Based on the stellar parameters for II Peg in Berdyugina et al. (1998) and us- 
ing the bolometric corrections in Flower (1996), we estimate the bolometric luminosity, the 
power emitted by the two stars over all wavelengths, to be Lboi ~5.5xl0^^ erg s~^. 

II Peg was first detected as a fiaring X-ray source by the Ariel- V Sky Survey Instrument 
(1.5-20 keV) in the mid-1970s (Schwartz et al. 1981) which detected 2 fiares with peak fiuxes 
of 1.1 X 10^^ and 1.5 x 10~^ erg cm^^ s^^ in the 2-10 keV band, equivalent to peak X-ray 
luminosities of 1.6 x 10^^ and 2.8 x 10^^ erg s'^. Since then, moderate and large flares from II 
Peg have been detected by many subsequent X-ray and EUV observatories, including Ginga 
(Doyle et al. 1991), EXOSAT (Taghaferri et al. 1991), EUVE (Patterer et al. 1993; Osten 
& Brown 1999), ASCA (Mewe et al. 1997), BeppoSAX (Covino et al. 2000), and Chandra 
(Huenemoerder et al. 2001). 

In this paper, we report on a large flare observed on II Peg by the detectors on the 



This is in contrast with the solar corona, where the coronal iron abundance appears to be enhanced by 
a factor of ^4 over the solar photospheric abundance (Feldman & Laming 2000). For more discussion on 
stellar coronal abundances, see Giidel (2004). 
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Swift Gamma- Ray Burst Mission (Gehrels et al. 2004). Following the discussion in Schaefer 
et al. (2000) and the characteristics of the flare (radiated energy, luminosity) we term this 
event a "superflare". The event triggered Swift's Burst Alert Telescope (BAT) in the hard 
X-ray band on December 16, 2005 at 11:21:52 UT^, and its characteristics in the soft and 
hard X-ray bands were determined for the extent of time for which hard X-ray emission 
was detected, ~ 7000 s. The UV/Optical Telescope did turn on during the trigger, but the 
images were saturated due to II Peg's brightness (Vf5i7.4), and the data were of little use. 
In this paper we concentrate on analysis of the X-ray Telescope (XRT) and BAT spectra. 



2. Data Reduction 

2.1. BAT Data Reduction 

BAT is a coded-mask instrument with a very large (~1.2 sr) field of view. The con- 
struction and operation of the instrument is discussed in Barthelmy et al. (2005). An 
overview may be found at: http:/ /swift.gsfc.nasa.gov/docs/swift/about_swift/bat_desc.html. 
All Swift data are made public as soon as possible, and the raw BAT data are available by 
navigating from the above site. 

The BAT data used in this analysis are "survey" data, intended for the BAT all-sky 
hard X-ray survey. It consists of accumulated counts, in time bins from 60 to 300 seconds, 
for each of the 32k CZT detector elements, in 80 energy bins nominally from 10 to 194 keV 
plus high and low integral energy bins. The survey data are supplemented during an interval 
around the BAT trigger time, including during the Swift slew maneuver, by "event" data, 
which contain higher resolution time and energy data on each photon count in the detector 
array. 

We used BAT pipeline software within FTOOLS^ version 6.0.3 to correct the energy 
from the efficient but slightly non-linear energy assignment made on board. For the spectral 
data reported here, we used batbinevt to produce mask-weighted spectra in several broad 
time intervals. For the light curve data we created sky images in two broad energy bins for 
each time interval using batbinevt and batfftimage, and found the fiux at the source position 
using batcelldetect, after removing a fit to the diffuse background and the contribution of 



^http: / /gen. gsfc.nasa.gov /gcn/gcn3 /4357.gcn3 

^The FTOOLS software paekage provides mission-specific data analysis procedures; a full description of 

the procedures mentioned here can be found at http://heasarc.gsfc.nasa.gov/docs/software/llieasoft/; the 
names of individual procedures referred to in this and the following section are italicized. 
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bright sources in the field of view. 

As seen in Figure 1 , the spectral data are reported in 3 intervals relative to the BAT 
trigger time of 11:21:52 UT: "Trigger", T-760 to T+126 sec; "Orbit 1", T+126 to T+1661 
sec; and "Orbit 2", T+4810 to T+7451 sec. There was no significant fiux in the BAT energy 
range during Orbit 3. To increase the signal-to-noise ratio at energies above 38 keV, the 
BAT spectra were grouped by 3 bins. 

2.2. XRT Data Reduction 

For a technical description of XRT and its operations, we refer the reader to Burrows 
et al. (2005). XRT started observing at 11:24:16 UT, that is, 144 s after the BAT trigger. 
XRT observed IlPeg in three different snapshots in three following orbits. The first snapshot 
("Orbit 1") lasted 1518 s. The second snapshot ("Orbit 2") started at 12:43:28 and lasted 
2545 s. The last snapshot ("Orbit 3") began at 14:24:16 UT and lasted 535 s. To produce the 
cleaned and calibrated event files, the data were reduced by means of the xrtpipeline task and 
calibration files of the CALDB 20051221 release^. Since the source count rate was over 60 
counts per second during the entire observation all data were collected in windowed timing 
(WT) mode. In order to extract the spectrum and light curve, events were selected with grade 
0-2^ from a 20 pixel (47") box corresponding to 80% of the encircled energy fraction (EEF) 
at 1.5 keV. We restricted our analysis to the energy band 0.8-10 keV, ignoring lower energy 
channels. In fact, the effective area calibration files (ancillary response functions, or ARFs) 
included in the CALDB 20051221 release were still preliminary and present systematics in 
the energy range 0.2-0.8 keV, which are noticeable in this high intensity event. Background 
extraction regions were chosen to be the same shape as the source extraction region, at 
approximately 50 pixels from the center of the source. The mean count rate in the source 
region in the three orbits was 161, 118, and 79 counts per second respectively, whereas the 
background count rate never exceeded 0.5 counts per second. In order to study the spectral 
variations of the Fe K fine (see below), we split the first orbit in three segments (144- 
650, 650-1156,1156-1662 s from the trigger), and the second in two (4902-6174,6174-7447 s 
from the trigger). The light curve was produced using the standard task Icurve and with 
a temporal bin of 30 sec. The spectra were grouped with the task grppha in order to have 



description and history of calibration files for Swift can be found at 
http: / /swift .gsfc.nasa.gov /docs /licasarc /caldb / swift / 

^See discussion in §2.3 of http://swift.gsfc. nasa.gov/docs/swift/analysis/xrt^wguide_vl_2.pdf for classi- 
fication of events and grades. 
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20 counts as a minimum for each energy bin. The ancillary response function was produced 
by means of the xrtmkarf taking into account the point-spread function correction (Moretti 
et al. 2005). 

3. Spectral Analysis 

Figure 1 shows the 0.8-10 keV XRT light curve and BAT light curves in two energy 
bands. It is evident that the harder X-ray energies have progressively earher peaks: at 0.8- 
10 keV, the light curve peaks ~1480 s following the BAT trigger, while, for the 14-40 keV 
band, the peak is ~840s post trigger, and for the 40-100 keV band, the maximum count 
rate occurs only ~480 s after the trigger. This points to the hardest X-ray photons having a 
more impulsive behavior than the soft X-ray photons. To examine the origin and evolution 
of these hard X-rays, we examined spectra corresponding to four different time intervals: 
the trigger spectrum and data collected in spacecraft orbits 1, 2 and 3. Spectral fitting 
was performed using XSPEC (vl2.0)^. To describe the hue and continuum emission from 
a thermal, diffuse, coUisionally ionized plasma, we used a customized version of the APEC 
model (Smith et al. 2001)^ which has continuum emission calculated out to photon energies 
of 100 keV. Spectral fitting proceeds via a forward method, computing model spectra with 
varying input parameters, convolving these with the instrument resolution and sensitivity 
and comparing them with observed spectra until suitable statistical convergence is achieved. 
We allow elemental metal abundances to vary in fitting each XRT spectrum, but because of 
the low spectral resolution in the soft X-ray region, the abundances are scaled to a common 
multiple A of the solar abundance (e.g. for a best-fit scaled abundance A =0.5, the Fe/H is 
2xFe/'lisoiar: ^-nd the Mg/H is also 2xMg/}isoiar)- We use the solar photospheric composition 
of Grevesse & Sauval (1998) in our spectral fitting; all fitted abundances quoted in this paper 
arc with respect to these values. In particular, the Fe/H ratio in Grevesse & Sauval (1998) 
is 3.16x10"^ 

There are deviations of the model from the data at low energies; the most prominent of 
these is a feature near 1.09 keV which appears in all XRT spectra. The residual near 1.09 
keV looks very similar, both in amplitude and shape, to those of the Crab spectrum used to 
calibrate the XRT. This feature is most likely an instrumental effect similar to other features 
observed in the (excluded) energy range below 0.8 keV. There are other large residuals around 



^XSPEC, an X-ray spectral fitting package, is part of the Xanadu software package, which is released 
packaged with the FTOOLS described previously and found at the same website. 

^For more information see http:/ /cxc. harvard.edu/atomdb/ 
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I keV; these may indicate that, instead of the Ne abundance scahng with the Fe abundance 
(in which case, the Ne/Fe abundance ratio would be one with respect to the solar ratio, 
which is 3.8), the Ne/Fe abundance ratio is greater than this value. Alternatively, it could 
also indicate the presence of a lower temperature component, which we could not constrain, 
as its major contributor is at energies <1 keV. 

Wc have also included photoelectric absorption in the spectral fit, with the column 
density Nh as a free parameter. The best-fit values do not show significant variation, but 
are systematically higher than previously determined interstellar column densities towards 

II Peg ('-^SxlO^^ cm~^; Mewc ct al. 1997; Patterer et al. 1993): this discrepancy is likely due 
to calibration uncertainties in the XRT at low energies. 

3.1. Trigger Spectrum 

The BAT trigger spectrum could be fit with a single component model, either a thermal 
(APEC) model with abundances fixed to solar (Grevesse & Sauval 1998), or a power-law 
model. At the high temperatures returned by the APEC model (~180 xlO^K) continuum 
emission dominates. Both models fit the data statistically equally well; however, we favor 
the thermal model based on previous behavior in large stellar flares (Favata & Schmitt 
1999; Pallavicini 2001; Franciosini et al. 2001) which behavior indicate the presence of high- 
temperature thermal plasma. Spectral flt results are given in Table 1 and shown in Figure 2. 

3.2. Orbits 1 and 2 

We investigated the presence of large residuals around 6.4 keV in the XRT spectrum 
during Orbits 1 and 2, and added a fine at this energy, flxing the energy at 6.4 keV and 
width at the instrumental value. The right bank of panels in Figures 3 and 4 display the 
5.5-9 keV region in Orbit 1 and Orbit 2 which includes this feature. We interpret this excess 
to be the signature of Ka emission from neutral or low ionization states of iron, and we defer 
the discussion of its origin and interpretation to §4.5. 

Only in two time intervals (Orbits 1 and 2) were both XRT and BAT spectra available, 
covering the 0.8-200 keV energy range with significant signal to noise. Our spectral fltting 
for these times started with a single thermal component, adding in additional thermal com- 
ponents to reduce the fit statistic. The large value of the second temperature component 
can be confirmed by examination of the 5.5-9 keV region (right bank of panels in Figures 3 
and 4) which contains several line diagnostics sensitive to 40-160 xlO^ K plasma. The 
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He-like and H-like Iron transitions at ~6.7 and ~6.9 keV are the most prominent, and their 
ratio constrains there to be plasma at ~ lO^K. Since the single temperature fit was not 
adequate to the spectra in either Orbit 1 or 2, we do not discuss it further. Temperatures 
near 10^ K have also been inferred from previous large stellar flares (e.g., Maggio et al. 2000; 
Franciosini et al. 2001). Figure 3 shows the spectral fltting results in Orbit 1, and Figure 4 
shows the spectral fitting results for Orbit 2. The left bank of panels displays the fits to 
the 0.8-200 keV energy range for two temperature (2T+G), three temperature (3T+G), and 
two temperature plus nonthermal (2T+G+NT) models (all have a Gaussian component at 
6.4 keV as described above). The residuals plotted are deviations of the model from the 
data, calculated in units of a. The two-temperature fit has several energy bins above 40 keV 
(Orbit 1) and 20 keV (Orbit 2) with significant residuals, which prompted us to add in the 
additional component. 

In the case of three thermal components, the highest temperature component is mod- 
elled using a thermal brcmsstrahlung model whose contribution in the spectrum appears 
mainly at high energies. The normalization of the bremsstrahlung component is given in 
neUidV, as compared to the APEC model normalization neUndV, rii being the ion number 
density and Hh the hydrogen number density. For both orbits, the xl statistic is mini- 
mized with three thermal components, but the temperature of the highest component, ~ 
300 xlO^K, in both Orbit 1 and 2, is very high, and we reject this model for the following 
reasons. At high temperatures, the conductive cooling time dominates. We can compute 
the ratio of the thermal relaxation time of the plasma to the timescale for conductive cooling: 

Trelax/Tcond = 200-2-y2 (l) 

^10-^9 

(see discussion in Benz 2002), where Tg is the temperature in units of 10^ K, nio is the 
electron density in units of 10^° cm~^, and Lg is the loop length in units of 10^ cm. For 
T8=3, the timescale over which the plasma can relax to a thermal distribution at this tem- 
perature exceeds by a large factor the timescale on which the plasma would lose its energy 
via conductive losses, unless the density is very high and/or the length scales involved are 
very long. At high densities, the magnetic field required to confine the plasma becomes large, 
Bcon/=V^87m^A5^T = GOy/riioTsG. At low densities, the length scales exceed the primary stel- 
lar radius (~3R0) and the binary separation (asini =3.4x lO^'^^cm, Strassmeier et al. 1993). 
There have been numerous discussions in the literature concerning the thermal or nonther- 
mal nature of hard X-ray emission from solar fiares (Vilmer 1987; Brown & Smith 1980). 
Because other fines of evidence exist which confirm that particle acceleration does occur dur- 
ing stellar fiares (mostly microwave gyrosynchrotron emission), we consider a model which 
describes bremsstrahlung emission from suprathermal electrons. The arguments which apply 
to 3x10^ thermal plasma apply also to the second temperature component found in Orbits 
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1 and 2 (T2 from 1.2-1.5 xlO^ in Orbit 1, T2 from 7-9 xlO^ K in Orbit 2) with only a slight 
easement of the physical restrictions described above. However, the line information present 
in the He-like and H-like iron transitions confirm the existence of this plasma, whereas the 
3x10® component comes from continuum emission at hard X-ray energies. 

Thus, we consider an additional model which contains a nonthermal component. This is 
based on hard X-ray emission from solar flares, which usually shows evidence of bremsstrahlung 
radiation from a population of nonthermal electrons (Dennis & Schwartz 1989). If the emis- 
sion arises as the result of a beam of injected electrons propagating downward through the 
atmosphere, we expect the observed radiation spectrum to differ from the injection spec- 
trum, being modified by collisions with the increasingly dense atmosphere. The formulation 
for describing such "thick-target" hard X-ray emission is described by Brown (1971), where 
the observed spectrum 7(e) has the form 

/(e) ~ e''^ photons cmT^ keV~^ (2) 

where e is the photon energy and 7 is the photon index. Under the formalism of thick-target 
bremsstrahlung emission, the observed spectrum can be related to the spectrum of injected 

electrons, F{E) ~ erg s~^ keV~^, where E is the electron energy, and S is the power-law 
index of the electron distribution function; d is related to 7 as 5 = 7 -|- 1. We implemented 
a thick-target bremsstrahlung code® in XSPEC to model the behavior of the hard X-ray 
photons. The spectral shape is determined by the power-law index of the electron injection 
spectrum; the normalization depends on the total power input in accelerated electrons and 
the low-energy cutoff in the electron injection spectrum. We cannot constrain the cutoff 
with our dataset, and fix this value to an arbitrary value of 20 keV often used in solar flare 
analysis (Dennis & Schwartz 1989). Recent RHESSI solar flare results (Holman et al. 2003) 
have determined 37 keV as the highest value consistent with their data, although much lower 
energy cutoffs (down to a few keV) have also been determined (Sui et al. 2006). This value 
(20 keV) is roughly comparable with the approximate location in the spectrum where the 
nonthermal spectrum starts to dominate. Note that the choice of low energy cutoff has 
dramatic consequences for the power input in accelerated electrons: decreasing the cutoff to 
a lower value Ec increases the power input necessary to reproduce the spectrum by a factor 
{20/E^Y-^, or a factor ^ two higher for 5 ~ 3 and E^ = 10 keV. 

The fact that the 3rd temperature component in Orbit 2 is also unphysically high, 
coupled with the large residuals above 20 keV between model and data for a two-temperature 



^Based on the BREMTHICK code developed for analysis of RHESSI solar flare data available at 
http: / /hesperia.gsfc.nasa.gov/hessi/modelware.htm 
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fit and the xt statistic for the 2 temperature plus nonthermal model being nearly the same as 
for the three temperature model, leads us to conclude that there is evidence for nonthermal 
emission into the decay phase of the flare. Indeed, since the thermal plasma has cooled (T2 ~ 
8x lO^K in Orbit 2 compared with T2 ~1.4x lO^K in Orbit 1), the discrepancy between model 
and data for the two-temperature model shows up at smaller X-ray energies (see Figure 4). 
Based on either the three temperature model or the two temperature plus nonthermal model, 
in Orbit 1 the 0.8-10 keV flux is 6.12 xlO^^ erg cm'^ s'^ [(6.09-6.13) xlO-^ la uncertainty], 
and the 10-200 keV flux is 3.9 x 10"^ erg cm'^ g-^ [(3.6-4.0) x 10'^ la uncertainty]. Expressed 
as a fraction of the bolometric luminosity (Lboi ~5.5xl0^^ erg s~^), these values are 0.24 
(0.8-10 keV) and 0.15 (10-200 keV). In Orbit 2, the 0.8-10 keV flux is 3.62x10"^ erg 
cm-2 s-i [(3.61-3.64) X 10-9 la uncertainty], and 10-200 keV flux is 1.29 xlO"'^ erg cm'^ 
s~^ [(l.ll-1.37)xl0~9 la uncertainty]. Expressed as a fraction of the bolometric luminosity, 
these values are 0.14 (0.8-10 keV) and 0.05 (10-200 keV). 



3.3. Orbit 3 

In the third orbit following the trigger, there was no BAT detection of the source. The 
XRT did collect data for an interval of ~ 535 s, which we analyzed. There is no statistical 
evidence indicating the need for either a power-law component at the highest energies, or for 
a 6.4 keV feature. The thermal component is still at a relatively high level, and the 0.8-10 
keV flux is ~40 times higher than recorded during quiescent intervals seen with Chandra 
(5.2xl0~^^ erg cm~^ s~^, 0.5-12 keV; Huenemoerder et al. 2001). We conclude that the flare 
in soft X-rays lasted longer than the ~12 ks in which it was observed by Swift. Spectral flt 
results are listed in Table 1 and the data and model are plotted in Figure 2. 



4. Discussion 

The soft X-ray flux peaked during Orbit 1, with the radiated energy flux > 100 times that 
observed during previous observations made when the system was quiescent (Huenemoerder 
et al. 2001). Thus, we apply the "superflare" moniker to describe this event. There have been 
previous detections of stellar flares in the hard X-ray band of 20-50 keV (Favata & Schmitt 
1999; Pallavicini 2001; Franciosini et al. 2001; Krivonos et al. 2005), yet none have shown 
unambiguous evidence for nonthermal emission. Indeed, detailed spectral fitting for these 
fiares has indicated the presence of a supcrhot thermal plasma component (T>80 xlO^K), 
but, with no detections above ~50 kcV there has been no evidence of a power-law tail in the 
spectrum due to nonthermal emission. These results were somewhat puzzling, in that radio 
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observations necessitate continuous acceleration of electrons to describe the steady levels of 
microwave emission (Chiuderi Drago & Franciosini 1993). In addition, intense radio flares 
reaching 2-3 orders of magnitude enhancement over the quiescent level of microwave emission 
have been observed (Richards et al. 2003) which probably correspond to the "superflare" 
classification as well. 

We estimate the peak flux during Orbit 1 as if this were a solar fiare - the 1-8 A(1.55- 
12.40 keV) flux at Earth would be 440 W m^^ for d = lAU instead of d = 42pc. Using 
the notation for X-class solar X-ray flares {Xn = nxlO^'^ Wm^^), this would correspond to 
X4.4xlO*^. For comparison, the largest solar flare yet observed has been ~X30. Given the 
huge difference in radiative output, we should not expect this ffare to look and behave like 
solar flares. The key features of this superflare are (1) its intensity in the soft X-ray band, 
notably the evidence for high temperature thermal emission; (2) evidence for nonthermal 
emission at hard X-ray energies; and (3) the presence of 6.4 keV Iron K emission. We 
discuss the imphcations of each of these for the physics of flares in the following sections. 



4.1. High Temperature Thermal Emission and the Neupert Effect 

High temperature plasma (> SxlO^K) is evident in all four time intervals. The hottest 
component, l.SxlO^K, appears during the trigger, although this result is based only on 
10-200 keV continuum emission. The spectral fitting to Orbits 1 and 2 involving models with 
two or more temperature components returns a high temperature (~70-150 x lO^K) which is 
generally hotter in Orbit 1 than in Orbit 2. We have argued in a previous section (§3.2) that 
the temperature returned by modelling the spectra in Orbits 1 and 2 with three temperature 
components is not as plausible as the existence of a nonthermal thick-target bremsstrahlung 
emission component at hard X-ray energies, and so we will ignore the three temperature 
fits in favor of the nonthermal interpretation. The two temperature fit to the spectrum in 
Orbit 3 also reveals evidence for plasma at SxlO^K. Thus, the general trend appears to 
be that the plasma is cooling during this progression from the trigger through Orbit 3, as 
deduced by the decrease in the hottest temperature component from these spectral fits, with 
an exponential decay time of ~3 hours. 

The correspondence in many solar flares between thermal coronal energies and the 
amount and timescale of nonthermal energy deposition (the "Neupert Effect" ; Neupert 1968) 
lends credence to the posited physical association between particle acceleration and coronal 
heating. Stellar flares also sometimes show evidence of the Neupert effect, although typically 
radio emission or other proxies for hard X-ray emission (e.g. optical white-light ffares) are 
used to compare with soft X-ray thermal radiation (Giidel et al. 1996; Giidel et al. 2002; 
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Hawley et al. 1995). We have investigated the temporal relationship between the soft X-ray 
emission and hard X-ray emission in this flare to see if it shows the Neupert effect. In this 
case, the time rate of change of the thermal energy content of the corona (as revealed by 
soft X-ray radiation) should be roughly equivalent to the instantaneous nonthermal energy 
luminosity, Eth oc Lnt- The 40-101 kcV energy band of the BAT spectrum appears to be 
dominated by power-law hard X-ray emission, so we use this light curve to describe the 
temporal variation in nonthermal energy deposition. The thermal coronal energy is released 
mainly in the 0.8-10 keV energy band, so we use this light curve to constrain the tem- 
poral evolution in the thermal energy. Figure 5 displays the correspondence between the 
derivative of the XRT light curve and the instantaneous nonthermal hard X-ray emission. 
The greatest correlation occurs during the first part of Orbit 1. We take this as suggestive 
of a Neupert cfi'cct relationship between the nonthermal hard X-ray emission and thermal 
coronal radiation, although there is complex behavior here, as the nonthermal hard X-ray 
emission appears to persist into the decay phase of the flare. An additional complication lies 
in whether the plasma heating to such high temperatures takes place concurrently with, or 
as a result of, the particle acceleration producing the nonthermal emission. 

The Neupert effect argues for a causal relationship between particle acceleration and 
plasma heating. The energy dissipation which sets off the chain of events seen in a flare can 
heat plasma as well as accelerate particles (Dennis Sz Schwartz 1989); the exact proportion 
depends on the partition between plasma heating and particle acceleration, which is un- 
known. It is possible that plasma heating to the high temperatures seen during the trigger 
and subsequent orbits occurs before or concurrently with particle acceleration, although we 
cannot constrain the presence of the nonthermal component during the time of the trig- 
ger, due to poor signal-to-noise constraints, and we also cannot deduce what the pre-trigger 
conditions were. Bcnka & Holman (1994) discussed a model in which both thermal and 
nonthermal emissions occur and arc physically linked due to electron currents both heating 
plasma by Joule dissipation and accelerating electrons via a runaway process. Secondary 
heating such as predicted by the Neupert effect is then a different effect. McTiernan et al. 
(1999) show that high temperature solar flare plasma {T > 16xlO^K) is more hkely to ex- 
hibit the Neupert effect than low temperature plasma, but this assumes that the flare energy 
release occurs predominantly in nonthermal electrons. Li et al. (1993) computed hard and 
soft X-ray time profiles using models in which the hard X-ray emission was produced by 
either a super-hot thermal component or nonthermal thick-target bremsstrahlung, and con- 
cluded that the thermal model fails to reproduce the derivativity relationship between hard 
X-ray and soft X-ray emission. A situation where primary plasma heating and secondary 
plasma heating due to energy lost by nonthermal electrons would be more complicated. The 
current data appear to be consistent with a Neupert effect relationship, but a more confident 
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interpretation relies on knowing the times when nonthermal emisison and high temperature 
flare emission appeared. 

As discussed by Feldman et al. (1996), there appears to be a relationship between 
flare temperature and emission measure for solar flares and large stellar flares; "bigger" 
flares (more intense, larger emission measure) tend to be hotter. The high temperature 
component of the flare under consideration here also appears to fit with this general relation. 
Battaglia et al. (2005) found a correlation between nonthermal hard X-ray flux and thermal 
plasma parameters from a sample of solar flares, indicating that flares with large values of 
nonthermal emission also have higher temperatures and emission measures. Applying this 
result to the current flare confirms that a strong bias exists in detecting nonthermal hard 
X-ray emission from stellar flares, as very large and energetic (and hence rare, due to the 
flare frequency-energy relationship) flares are needed to achieve detections of nonthermal 
emission at hard X-ray energies. This problem is further compounded by the addition of 
the superhot thermal component, itself a consequence of the large flare, which complicates 
detection of the nonthermal component in the hard X-ray spectrum. 



The power Fq in the accelerated electrons as deduced by thick-target bremsstrahlung 
spectral flts is ~ 10^^ erg s~^, as hsted in Table 1. As noted above, this number depends on the 
value of the low energy cutoff in the injected electron spectrum, which we cannot constrain. 
However, our somewhat arbitrary choice of 20 keV as a lower energy cutoff does fit with the 
region of the spectrum where nonthermal emission begins to dominate. The spectral indices 
of the electron distribution returned from the spectral fitting are 6 of 2.8 and 3.1; these 
compare favorably with spectral indices inferred from solar HXR burst spectral indices (e.g., 
McTiernan & Petrosian 1991); significant evolution of the hard X-ray spectrum is usually 
observed during solar flares, but our sensitivity constraints do not permit examinations on 
a flner time scale. 

The flare lasted for > 11000 s in the soft X-ray band, but the timescale for hard X-ray 
emission is shorter, ~ 7000 s, based on the nondetection by BAT in the third orbit follow- 
ing the trigger. The major energy loss mechanism for the nonthermal electrons will be via 
collisions with the ambient thermal electrons; we thus estimate the energy loss timescale as 



where EkeV is the energy of the electron in keV, Ue is the ambient plasma density in the 



4.2. Nonthermal Hard X-ray Emission 




(3) 



-14- 



region of energy loss in cm~^, and In A is the Coulomb logarithm, falO under typical chromo- 
spheric conditions. We use rig of ~10^°-10^^ cm~^ based on electron density measurements of 
II Peg's lower transition region/upper chromosphere given in Doyle et al. (1992). Under these 
conditions, a 100 keV electron will lose its energy in 2-20 s. This is significantly shorter than 
the observed timescale for the hard X-ray emission, and thus the observations require either 
continuous acceleration and/or the presence of much more energetic electrons than can be 
diagnosed with the spectral energy coverage and sensitivity of the Swift detectors. The first 
option is broadly consistent with hydrodynamic stellar fiare models which require continued 
heating (Reale et al. 1997), hence continuous particle acceleration, assuming the two pro- 
cesses are physically finked. It is likely that both continuous acceleration and a population 
of highly energetic electrons exist; radio observations of active stars indicate the presence of 
~ MeV electrons. We expect that had radio observations of this flare been obtained, the 
observations would have shown a significant enhancement above typical levels. Indeed, using 
the Lx-^iR relationship of Giidel & Benz (1993), we would expect a centimeter- wavelength 
radio luminosity of (4.5-45) xlO^^ erg s^^ Hz"\ or 0.2-2 Jy flux density, using the 0.8-10 
keV X-ray luminosity during Orbit 1 (Lx=7.7xl0^^ erg s"^). Although, in most solar flares, 
particle acceleration occurs only during the impulsive phase, there is evidence for continued 
particle acceleration and chromospheric evaporation during the gradual phase of some flares 
(Cfiver et al. 1986; Kai et al. 1986; Czaykowska et al. 1999), which are typically long-duration 
events. 



4.3. Occurrence Rate 

The Swift mission had been operating for roughly 9 months before this event was de- 
tected, observing ~60% of the sky each day. This suggests a rough frequency of occurrence 
for these superflares of once every 5.4 months, or roughly once in 164 days. Other flare events 
from active binary systems with high luminosities have been only sporadically reported in 
the literature. If such events follow the usual power-law dependence of flare frequency with 
luminosity, then this is almost certainly due to the low probability of catching such a flare 
during pointed observations lasting a few days or less. For the 66 days of exposure accu- 
mulated in the >lneZ- F observations of Schwartz et al. (1981), the rate of flares in excess of 
6 X 10~^° erg cm~^ s~^ ( 10 — 30 times the typical 'quiescent' X-ray flux of 2 — 6 x 10~^^ erg 
cm~^ s~^ since determined for this system) in the 1.5-20 keV band was inferred to be 11 
per year. 

Dedicated surveys obviously will have a higher probability of detecting these events. 
Waldram et al. (2003) have a Ryle Telescope 15 GHz light curve of II Peg (their Fig. 13) 
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which spans 14 months, albeit with some gaps, which shows one flare with an intensity 
greater than 200 mJy (fs 100 times the quiescent radio flux density) and 2 others >100 
mJy. Based on several years' observation at centimeter wavelengths with the Green Bank 
Interferometer reported in Richards et al. (2003), we can also estimate the occurrence rate of 
radio supcrflares in active binary systems, using the flux data for HR 1099 and UX Arietis ^, 
two analogous active binary systems to II Peg containing a K subgiant and hotter companion 
(but where the K subgiant, as in II Peg, is presumed to dominate in the X-ray and radio 
emission). Only for 0.2% of the time did flare events in HR 1099 and UX Ari reach radio 
luminosities in excess of 100 x the quiescent values, corresponding to a total of 9 flare events 
for these two systems. The average time between superflares is 136 days (165 days for 4 
events on HR 1099 and 113 days for 5 events on UX Ari); with the usual caveat of small 
number statistics, the event rate for X-ray superflares and radio superflares appear to be 
broadly consistent, being roughly once every several months to a year. 



4.4. Contribution of Thermal and Nonthermal Energies 

A lower limit to the amount of thermal energy in the plasma can be obtained by de- 
termining the radiative flux from the thermal plasma over a large range of photon energies. 
Wc did this in XSPEC through the use of a "dummy" response covering 0.01-200 keV. For 
Orbits 1 and 2, we used the best-fit two-temperature plus Gaussian plus nonthermal model 
parameters listed in Table 2, and removed the Gaussian and nonthermal model components. 
A lower limit to the thermal energy can then be calculated as iircPFAt, where F is the 
radiative flux from 0.01-200 keV (9.5x10^^ erg cm~^ s~^ for Orbit 1 and 4.8x10"^ erg cm~^ 
s"^ for Orbit 2), and At is the duration of Orbit 1 and 2 (1535 and 2641 s, respectively). 
In both cases, the radiative energy estimates work out to ^ 3x10^^ erg, for a total of 
6x10^^ erg. The total radiated energy of the hot plasma over all wavelengths can also be 
computed from the radiative loss function, using the radiative losses at the temperatures 
returned from spectral fitting appropriate for the derived abundances, and multiplied by the 
emission measure at that temperature and total duration. Wc used the line and continuum 
emissivities in APEC to calculate this function. The total radiated energy calculated in this 
manner matches to within an order of magnitude the energy estimates above; this is due 
partly to the wide wavelength range considered in the "dummy" response. We note that 
the total radiated energy over all wavelengths from the few solar flares in which changes in 
total solar irradiance have been made (Woods et al. 2004) indicate that the total radiated 
flare energy can be larger than the radiated energy of the hot plasma by a factor of up to 



^Obtained from ftp: / /ftp. gb. nrao.edu/pub/fghigo/gbidata/gdata/gindex. html 
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ten. Thus, a lower limit to the radiated energy for the flare in Orbits 1 and 2 combined can 
be placed at ~ 10^^ erg. This is the amount of energy the plasma loses by radiation, and 
ignores the effect of conductive energy losses (which will be important at high temperatures) 
as well as bulk kinetic energy of the plasma or any energy loss due to expansion. 

We can estimate the amount of conductive energy losses by using the equation for con- 
ductive flux, 

i^cond = «7^'/'VT = — ^ erg cwT'' s"^ (4) 

where k is the Spitzer conductivity value (= 8.8 x 10"''' erg cm"-*^ s~^ K"'''/^), and we have 
approximated VT" ~ T/L where L is a characteristic length along which conductive energy 
is lost. The energy density per unit time can then be estimated as F^^^jL ~ kT'^I'^ jl?. 
Using the relationship between volume, volume emission measure, and electron density 
{VEM = nIV), and the duration of each orbit, an estimate of conductive energy losses 
in Orbit 1 and 2 can be estimated as 

KT^^VEMAt 

-Econd = ^'^9 ■ (5) 

Using the temperature and volume emission measure^*^ the conductive energy lost in Orbit 

1 is ~ 4xlO^^/(L9nfo) erg, while in Orbit 2 it is ~6xlO^^/(L9nfo) erg, for a length in units 
of 10^ cm and electron density in units of 10^° cm~^. This is probably an upper limit on the 
amount of conducted energy, as evidence from solar flares indicates (Jiang et al. 2006) that 
the classical expression for conductive flux may overestimate the total energy decay rate. 

We estimate the total energy represented by nonthermal electrons during orbits 1 and 

2 by multiplying the power output by the exposure time. We cannot diagnose variations 
in the acceleration of electrons on shorter timescales, due to signal-to-noise limitations, so 
these calculations assume that electron acceleration occurs at a constant level during each 
orbit. In Orbit 1 Fq x At =10"^° erg, and in Orbit 2 Fq x At =2 xlO'^^ erg, where the low 
energy cutoff of 20 keV was used. In order to match the thermal energy (here assumed to 
be dominated by conductive energy losses) with the nonthermal energy input, we obtain a 
constrain on the product of electron density and length scale, LgUio =60 for Orbit 1 and 
LgUio =20 for Orbit 2. We have no independent constraints on either quantity. We note that 
the reasonable agreement between the nonthermal energy and the upper and lower limits on 
the thermal energy estimates add support to the nonthermal interpretation. 



J HeUHdV tabulated in Table 2 can be related to the volume emission measure n^V by using He/nn =1.2. 
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4.5. Iron Fluorescence 

There is excess emission redward of the Fe XXV and Fe XXVI features visible in the 
XRT spectra of Orbit 1 and 2, visible in the right-hand bank of panels in Figures 3 and 4. 
We attribute this to emission from the iron Ka feature at 6.4 keV. Both electron-impact 
ionization or photoionization mechanisms are capable of removing K shell electrons from 
neutral or near-neutral iron in the photosphere and producing the 6.4 keV spectral line 
(Emshe et al. 1986), although no solar Ka line has been definitively identified with electron- 
impact ionization (Parmar et al. 1984). We subdivided the spectra in orbits 1 and 2 to 
investigate the time evolution of the 6.4 keV feature in more detail. Figure 6 displays a 
close-up of the 5.5-8 keV region, and Table 3 lists the results of spectral fitting for the five 
different time-resolved spectra. 

Previous observations of Ka emission from stars other than the Sun have been of pre- 
main sequence stars and have concluded that the fluoresced material is located in the circum- 
stellar disk which is bathed in the thermal hard X-ray continuum radiation emitted by flaring 
plasma. Applying this interpretation to the event on II Peg poses several problems, however. 
The equivalent width of the 6708 Alithium resonance line and the lithium abundance are 
not consistent with a pre-main sequence evolutionary state, nor is the observed C/N ratio 
(see discussion in Berdyugina et al. 1998). The system is therefore generally considered to 
be an old-disk population star. On the other hand, an infrared excess appears to have been 
detected in II Peg (Lazaro et al. 1987; Rodono et al. 1998) which is rather unexpected for a 
Pop I star. Further, the column densities returned from X-ray spectral fitting are systemat- 
ically higher than previously determined interstellar values for this line of sight, suggesting 
the possibility of additional absorption from a circumstellar disk. However, we consider the 
column density returned from spectral fitting highly suspect, as we have excluded the en- 
ergy range below 0.8 keV due to systematics, and there could be additional unrecognized 
systematics affecting the absorption column density. Thus we consider it is more plausible 
that the 6.4 keV signature arises from photospheric, not circumstellar, material, although 
we cannot at the moment rule out a face-on circumstellar disk. 



4-5.1. Fluorescence 

If we attribute the formation mechanism to a photoionization fluorescence mechanism, 
then the continuum radiation above the iron K edge, 7.11 keV, is the source of the photoion- 
ization producing the 6.4 keV feature. The luminosity and equivalent width of the feature 
can be used to deduce the scale height of the coronal emission above the photosphere, using 
equation 2 in Tsujimoto et al. (2005) 
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where Ai7 is the sohd angle subtended by the photosphere as seen by the flaring X-ray 
source. The continuum spectrum in this region is dominated by the hot thermal plasma, 
with a spectral shape I{E) oc exp~^/*^^. We use the fluorescence yield Yxa of 0.342 
(Bambynek et al. 1972). Ej^a=6.4 keV and x is the iron K edge energy, 7.11 keV. The 
quantity upeiE) is the photoelectric cross-section of iron; we use upeiE) — 2x 10~^°(£'/x)~^ 
cm^ (GuUikson 2001). We rewrite the integral involving npeis) as J {AFe)'nHds, where Ape — 
npe/nH- A comparison of coronal iron abundances (here and in Huenemoerder et al. 2001) 
and photospheric iron abundances does reveal chemical fractionation occurring in II Peg's 
atmosphere (in a sense opposite to that seen in the Sun), yet we cannot constrain such spatial 
variations. The hydrogen column density in the chromosphere/photosphere is estimated 
using the column mass density in the atmospheric models of Short et al. (1998). Between the 
temperature minimum and 10^ K, the column mass ranges from ~0. 001-2 g cm~^, or 10^^- 
2x10^^ cm~^ with a mean molecular weight ji ~0.6. Since the photospheric iron abundance 
of II Peg is higher than the coronal abundance by factors of 2-4, and the hydrogen column 
density lower in the atmosphere exceeds that in the corona by orders of magnitude, we 
assume that the major contribution to this integral occurs in the photosphere and express 
the integral as Apcphot f nnds = Ape.photNH, where Apcphot -0.4 kpcsoiar = 1.26 xlO"^ 
using the revised iron abundance of Grevesse & Sauval (1998), and the above ranges of Nh- 
In the solar case, the main contribution also occurs in the photosphere, due to the increasing 
column mass as one proceeds from the corona to the photosphere; downward propagating 
photons with E>7 keV are optically thin to photoelectric absorption and Compton scattering 
(see discussion in Parmar et al. 1984), and the photospheric iron abundance being ~one 
fourth that of the solar coronal iron abundance does not affect where in the atmosphere the 
main contribution to fluorescence occurs 

The sohd angle AQ can be expressed in terms of a height using the fraction of photons 
intercepted by the star, (Bai 1979), 



An{h) = 27r 



(7) 



R^ + h 

where h is the scale height of the flaring X-ray source; we take to be the radius of the 
primary of the system, a; SR© (Berdyugina et al. 1998). Figure 7 displays the relationship 



^^Note that in the solar case, with an independent constraint on the height of the soft X-ray source from 
e.g. X-ray imaging telescopes, the flux of the Ka line can be used to deduce the photospheric iron abundance 
(Bai 1979). 



-19- 



between height and column density for the measured equivalent width and plasma temper- 
ature. There is an asympotic dependence on column density, so that the minimum column 
densities which can reproduce the equivalent widths arc of order 10^^ cm~^. This would 
indicate, based on the atmospheric modelling of Short ct al. (1998), that (1) the fluores- 
cence occurred deep in the atmosphere, near the temperature minimum region, and (2) the 
maximum flare scale height is 0.5R*. There appears to be a discrepant behavior in the first 
sub-segment of Orbit 2, where the equivalent width shows an anomalous value compared to 
Orbit 2b and Orbit Ic. There are problems with this interpretation, however, due to the 
fact that at Nh ~ I/ctt ~ 1-5 x lO^^cm"^, where ar is the Thomson cross section, photons 
start experiencing significant Compton scattering. In contrast with the solar photosphere, in 
which the temperature minimum region is reached with a lower column mass density (Ver- 
nazza et al. 1981), photons of energy 6.4 keV will also experience significant photoelectric 
absorption at such high column densities, decreasing the efficiency of producing fluorescence 
emission in such an environment. Thus we conclude that the fluorescence mechanism is 
probably not a valid interpretation for the formation of the 6.4 keV line. 

The iron Ka emission feature can have a great utility during stellar flares, if further 
observations con&m the fluorescence mechanism. With an independent constraint on the 
flaring loop height, say from hydrodynamic flare modelling relating the loop height to the 
flare evolution in the T-VEM plane (Reale et al. 1997), the equivalent width of the 6.4 keV 
feature can be used to constrain the photospheric Fe/H value. Line and continuum emission 
in the soft X-ray spectrum naturally constrain the coronal Fe/H value, allowing a simulta- 
neous measurement of both with a single data set. This would be advantageous to studying 
chemical fractionation in stellar atmospheres, particularly as photospheric abundances of 
active stars are notoriously difficult (due to e.g. fast rotation and/or binarity). 



4-5.2. Collisional Ionization 

An alternative explanation raised for the production of the 6.4 keV feature seen in solar 
flares is the collisional ionization of K shell electrons by a beam of nonthermal electrons 
(Emslic et al. 1986). We used equation 11 in Emslie et al. (1986) to estimate the 6.4 keV 
line flux (in photons cm^^ s^^) assuming the production mechanism is collisional ionization 
by a beam of nonthermal electrons 

^ = W^7^^7^ ^^^'-^^' iV*))--^g.(^)ciS , (8) 

where u is the fluorescence yield of iron, /3 is the branching ratio between Ka and K/3, =0.882 
(Bambynek et al. 1972), d is the distance to the object, K — 27re^A/^^ with e the electronic 
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charge in e.s.u., A the Coulomb logarithm, d the conversion from keV to erg, = 1.6x10"^, 
7 is the photon spectral index (related to the spectral index of the electron distribution by 
7 = 5 — 1). The value Eiow is the low-energy cutoff of the electron distribution, fixed to 
20 keV in our spectral fitting. The integral extends from the Ka edge of 7.11 keV to infin- 
ity; Qi{E) is the collisional ionization cross section, obtained from the theory of Arthurs & 
Moiseiwitsch (1958). Eq is the initial energy of the electron at injection, and E is its energy 
after collisional (thick-target) encounters. Specifically, Eq satisfies the equation 

El - 3KN*Eo = E^ (9) 

where N* is the column density in the Kcc-emitting region. The results of the thick- 
target bremsstrahlung spectral modelling during orbits 1 and 2 are used to estimate the 
6.4 keV photon flux, as there is not enough signal to subdivide the BAT spectra as was done 
above for the XRT spectra. For npe/nH we initially used the photospheric iron abundance, 
^Fe=l-26xl0-^ 

We calculate the 6.4 keV line flux relative to the underlying continuum flux (determined 
from the unfolded spectrum in XSPEC) to determine the equivalent width of the feature. 
Table 2 lists the derived values for the results from thick-target spectral modelling in Orbits 
1 and 2. The right panel of Figure 7 displays the results for the thick-target model in 
Orbits 1 and 2 as a function of N*. The observed equivalent widths can be reproduced 
for values of N* <10^'^ cm~^. This mechanism is effective at lower column densities than 
a fluorescence mechanism imphes, and thus would take place higher in the atmosphere. If 
the abundance fractionation which is known to occur between the photosphere and corona 
in II Peg is happening at these column depths, the value of upe/nH appropriate for the 
calculations must consequently be lower. The equivalent widths will then be lower, by a 
maximum of ~ 1/4 the values in the Table 2, corresponding to the maximum iron depletion 
in the corona as found by Huenemoerdcr ct al. (2001). We conclude that the collisional 
ionization mechanism is to be preferred over the fluorescence mechanism. 

4.6. Loop Heights 

If we assume that the flare emission originated from an ensemble of coronal loops with 
uniform cross section and roughly semi-circular mid-plane shape, then we can express the 
observed volume emission measure in terms of the loop height, density, number of loops, and 
loop cross section. Following Equation (5) of Huenemoerdcr et al. (2001), the loop height is 

f N\ "^/^ f a\ -2/3 / VEM \ f Ue \ -2/3 ^ ^ ^ 
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where the volume emission measure (VEM) is rewritten as the contribution from N flaring 
loops, each with height h and aspect ratio a, electron density Ug. Huenemoerder et al. 
(2001) deduced loop heights of 0.05 Rj, under quiescent conditions, for N = 100, a = 0.1, 
VEM = 7.9 X 10^^ cm^^ and rig ~ 10^^ cm^^. The peak emission measure determined here 
is ^100 times larger than that in Huenemoerder et al. (2001); applying the same analysis to 
this superfiare yields h = O.SR^,, for the same density, for = 100. If instead only one loop 
is involved, h could be as much as 1.3 R^. An electron density higher than the adopted value 
would result in more compact loops. The fluorescence analysis of the 6.4 keV feature also 
yields a constraint on the solid angle extended by the photosphere as seen by the continuum 
X-ray source, i.e. the scale height of the corona above the photosphere. The scale heights 
inferred in §4.5.1 are consistent with this simple scaling. 



5. Conclusions 

We have identified nonthermal emission as the most plausible mechanism to explain 
the hard X-ray emission seen during the rise and decay of a large stellar flare. This flare 
interestingly also displayed evidence for emission at 6.4 keV whose formation we attribute 
to the mechanism of collisional ionization. The increased sensitivity of the Swift BAT has 
enabled a detection at much higher energies than had been possible with previous hard X-ray 
telescopes. Nonthermal hard X-ray emission has enabled an investigation of the energetics 
of large stellar flares, without complication from optical depth effects and source inhomo- 
geneities. The characteristics of this flare — a rare, intense, transient event — point to the 
value of triggered observations to study such flares; targeted observations would have had 
a very low likelihood of observing such an event during the short timescales over which the 
hard X-ray emission is produced. Triggered observing modes necessarily miss the preflare 
emissions, as the flux must pass a threshold value to warrant a telescope slew. Still, serendip- 
itous science can be obtained, as this Swift-observed flare demonstrates. Current plans are 
to increase the trigger sensitivity threshhold of Swift by a factor of 4.5, thus enabling more 
such opportunities. 

Multi-wavelength observations of solar flares have revealed the dynamical response of 
the atmosphere to these sudden intense inputs of energy. Future observations of stellar flares 
like the one discussed here can benefit from the global array of telescopes which have been 
harnessed for gamma-ray burst studies. It will hopefully be possible to take advantage of 
optical telescopes to reveal the response of the lower atmosphere and radio telescopes to 
explore further the action of nonthermal particles. Observations of such events will allow for 
a comparison of particle acceleration processes in active stars and the Sun. 
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Table 1. Trigger and Orbit 3 Spectral Fit Results ^ 



Parameter 


Trigger 


Orbit 3 




T-760:T+1262 


T+10944:T+114792 


Nh (1022 cm-2) 




0.095 






(0.07-0.12) 


Ti (lO'^K) 




15 






(14-16) 


{JnenHdV)i (10^^ cm-^) 




15 






(10-21) 


T2 (lO^K) 


177 


61 




(112-294) 


(54-70) 


\j ngTiHClV )2[^'J cm ) 




oZ. 




(11.-39.) 


(30.-36.) 




1. 


0.16 




fixed 


(0.12-0.22) 


X' (dof) 


23.4 (32) 


521 (514) 


/(0.8-10 keV)'' (erg cm'^ s"^) 




2.09x10-*^ (0.08)" 






(2.08-2.12) X 10-9 


/(10.-200. keV)^ (erg cm'^ s-^) 


1.2xlO-V^^-^^^/ 






(0.8-1.3) X 10-9 




At{s) 


834 


535 



^Unless otherwise indicated, error ranges refer to 90% confidence inter- 
vals. 

^Start and stop times of each time interval, in seconds since the trigger 
time of 11:21:52 UT. For the Trigger spectrum, these times refer to the 
BAT, for Orbit 3 the time is the XRT. 

^Scaled abundance relative to solar, using solar abundances of Grevesse 
& Sauval (1998). 

^Error ranges correspond to 68% confidence intervals. 
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■^Numbcr in italics refers to flux converted to luminosity and expressed 
as a fraction of the bolometric luminosity. 
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Table 2. Models Fit to XRT & BAT Spectra in Orbits 1 and 2^ 



Paranieters 


1T+G2 


2T+G2 


3T+G2 


2T+G+NT2 


— Orbit 1: T+126:Th-166P — 


N^^ (10^2 cm-2) 


0.062 


9.8e-2 


0.1 


0.1 




(0.057-0.067) 


(0.09-0.1) 


(0.09-0.11) 


(0.09-0.11) 




0.33 


0.33 


0.38 


0.34 




(0.28-0.38) 


(0.27-0.38) 


(0.31-0.49) 


(0.29-0.39) 


Ti (lO^K) 


139 


16.5 


15.4 


15.3 




(134-143) 


(15.0-18.4) 


(14.4-17.2) 


(14.4-16.0) 


(/ neUHdV)! (10^^ cm-3) 


85 


7.45 


5.3 


5.8 




(84-86) 


(5.7-9.4) 


(3.8-7.0) 


(4.5-7.0) 


(10"^ photons cm"2 s"^) 


2.4 


2.5 


2.6 


2.7 




(1.8-3.1) 


(1.9-3.2) 


(1.9-3.2) 


(2-3.3) 


EW^ (eV) 


45 


48 


47 


51 


T2 (lO^K) 




152 


118 


139 






(146-157) 


(102-132) 


(133-144) 


{JnenHdV)2 (10^4 cm-^) 




83 


65 


83 






(82-84) 


(47-77) 


(83-84) 


T3 (lO^K) 






306 










(224-580) 




/ rienidV (10^^ cm-^) 






20 










(7-38) 












2.8 










(2.4-3.4) 


Fo^ (10^6 erg s"^) 








8.5 










(7.6-28) 


(dof) 


1065 (836) 


943 (834) 


918 (832) 


924 (832) 




— Orbit 2: T+4 


810:T+745P- 






N^^ (1022 cm-2) 


0.011 


0.066 


0.067 


0.066 




(0.0055-0.016) 


(0.057-0.075) 


(0.058-0.077) 


(0.058-0.075) 


A^ 


0.31 


0.26 


0.28 


0.27 




(0.29-0.34) 


(0.24-0.29) 


(0.25-0.32) 


(0.25-0.30) 


Ti (lO^K) 


67 


15.5 


15.3 


15.4 
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Table 2 — Continued 



Parameters 


IT+G^ 


2T+G2 


3T+G2 


2T+G+NT 




(66-69) 


(15.2-16.4) 


(15.0-15.8) 


(15.1-15.8) 


(/ nenHdV)i (10^^ cm'^) 


58 


13 


12 


12 




(57-59) 


(12-15) 


(10-14) 


(11-14) 


(lO"'^ photons cm~^ s~^) 


1.14 


0.9 


0.95 


0.99 




(0.84-1.6) 


(0.56-1.3) 


(0.6-1.32) 


(0.65-1.3) 


EW (eV) 


42 


33 


34 


34 


T2 (lO^K) 




86 


74 


83 






(83-90) 


(67-81) 


(80-86) 


(J nenHdV)2 [10^ cm 




51 


47 


51.4 






('50-521 


(41-50) 


f50 7-51 51 


T3 (lO^K) 






302 

(191-1000) 




/ rien/dV^ (10^^ cm-3) 






6.3 

(2.5-12) 




5^ 








3.1 

(2.2-4.1) 


Fo^ (10^6 erg s'^) 








8.6 

(6.0-26) 


(dof) 


1460 (814) 


952 (812) 


932 (810) 


931 (810) 



^Uncertainty ranges refer to 90% confidence intervals. 

^Key to model combinations:nT= n thermal components, G=Gaussian at 6.4 keV, 
NT=nonthermal thick target bremsstrahlung model. See text for details. 

^Start and stop times of each orbit, in seconds since the trigger time of 11:21:52 
UT. 

''Abundance of thermal plasma as a multiple of the solar photospheric metal abun- 
dance given in Grevesse & Sauval (1998). 



"Flux emitted in the 6.4 keV line. 
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^Equivalent width of the 6.4 keV hne. 

^Power-law index of accelerated electron spectrum; see text for details. 

^Power in accelerated electrons from thick target bremsstrahlung modelling above 
a cutoff energy of 20 keV. 



Table 3. 6.4 keV Feature^ 



Parameter 


Orbit la 


Orbit lb 


Orbit Ic 


Orbit 2a 


Orbit 2b 




T+144:T+650^ 


T+650:T+1156^ 


T+1156:T+16622 


T+4902:T+6174^ 


T+6174:T+7447^ 






n 1 


n 1 1 


075 


048 




I'D 07—0 1 

I U.U ( U. lU I 


I'D nsT— n 19"! 

I u.uo t u. iz 1 


(0 DQ— n 1 9\ 


^Ti nfi5— n r\Q9\ 


I U.UOO U.UUl I 


Aa 


n 47 


n 35 


n 96 


98 


9Q 




lU.OU U.UOJ 




1 7-0 35"! 


9^0 39 






15 1 


15 3 


1 5 8 


16 2 


15 4 




iZ. / D— lo.OD J 


n ^7 1 8 7"! 

^^lo.o / — lo. / ) 








( f T) n zjdV)^ (1 n^*^ rTn~2'l 




6 n 


1 1 1 


1 5 1 


10 4 






8 8 y'i 

^^o.o— 0. ( J 


(7 A 1 7 7"! 


n 9 1 8 c:^ 
^^IZ.o— lo.O J 


('8 7 1 9 8'\ 
(^o. / — IZ.o ) 


T2 (lO^K) 


151 


142 


140 


88 


86 




(138-173) 


(130-155) 


(126-169) 


(82-95) 


(81-91) 


(/ nenHdV)2{lO^^ cm-^) 


6.8 


8.7 


9.3 


5.3 


4.9 




(6.6-6.9) 


(8.5-8.8) 


(9.1-9.4) 


(5.1-5.5) 


(4.8-5.0) 


Flux in 6.4 keV feature 


2.5 


2.7 


2.8 


0.65 


1.68 


(10~^ photons cm~^ s~^) 


(1.4-3.6) 


(1.5-3.9) 


(1.5-4.0) 


(0.08-1.1) 


(1.2-2.2) 


EW (eV)^ 


55 


48 


47 


18 


61 


(dof) 


526 (552) 


549 (538) 


592 (555) 


595 (579) 


621 (555) 



^Error ranges refer to 90% confidence intervals. 

^Start and stop times in Orbits 1 and 2 during which XRT spectra were extracted, in seconds since the trigger time 
of 11:21:52 UT. 



^Abundance of thermal plasma as a fraction of the solar photospheric metal abundance of Grevesse & Sauval (1998). 
''Equivalent width of the 6.4 keV hne. 



^ 180 
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Time since trigger (s) 



Fig. 1. — Light curves from XRT, 0.8-10 keV (top) and two hard X-ray energy bands from 
the BAT (14-40, 40-101 keV; bottom). Time bins where SNR<1 are shown as downward 
arrows. Time is expressed referenced to the trigger of 11:21:52 UT on 16 December 2005. 
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Fig. 2. — Spectral fit to BAT trigger spectrum (left) and XRT spectrum in Orbit 3 (right). 
For the trigger spectrum a single temperature component is fitted, while for Orbit 3 two 
temperature components are fitted. Data with error bars are indicated by crosses. Histogram 
shows final model. For Orbit 3 the contribution of the lower temperature component is shown 
with a dashed line, while the contribution of the higher temperature component is shown 
with a dotted line. Residuals are differences between data and model in units of a. Spectral 
fits are tabulated in Table 1. 
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Fig. 3. — Spectral fits to XRT and BAT spectra in Orbit 1. The top panel of each sub-plot 
displays data and model, while the bottom panel plots the residuals between data and model. 
The residuals are in units of a and are plotted with unity error bars. In the upper panel 
of each sub-splot, crosses show spectral data points with errors. The histogram displays 
the final model with all components included. The contribution of each additive component 
is also shown: the first temperature component is delineated by a dashed line, the second 
by a dotted line, and the Gaussian at 6.4 keV by the thin solid line. The large residuals 
above 40 keV in the two temperature fit require the addition of another model component. 
In models with either a third temperature component or a nonthermal component, this 
additive component is indicated by a dashed-dotted line. The left bank of sub-plots shows 
the data from 0.8-200 keV, while the right bank of sub-plots zooms in on the Fe K region. 
The He- and H-like transitions of Iron confirm the 2nd temperature component in all cases. 
See Table 2 for fit parameters. 
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Fig. 4. — Same as Figure 3, but for Orbit 2. Note the residuals between data and model for 
the two temperature fit in the left-hand panel show a systematic excess above about 20 keV, 
which can be fit by either a third temperature or a nonthermal model. 
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Fig. 5. — Plot of 40-101 keV BAT light curve (heavy squares with error bars), along with 
scaled derivative of 0.8-10 keV XRT hght curve (filled squares connected by a hne) for the 
first 2500 seconds of the flare. The XRT hght curve has been rebinned to a coarser time 
sampling than original data. The correlation for roughly one thousand seconds after the 
trigger may be indicative of the Neupert effect relating energy deposition in nonthermal 
particles and subsequent plasma heating. 
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Fig. 6. — Intra-orbit variation of the spectral region around the Fe Ka 6.4 keV feature, 
for three time shces in Orbit 1 and two in Orbit 2. Time intervals are listed in Table 3 
Data points and errors are shown as crosses. The 6.4 keV component is drawn with a solid 
line, and the dotted line delineates the emission from the APEC model; the thick histogram 
describes the total model emission over this energy interval. The prominent features at 
6.7 and 6.9 keV are the He- and H-like transitions of iron, respectively. "Resid" refers to 
residuals between model and data in each energy bin, in units of a. See Table 3 for spectral 
fit results. 
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Fig. 7. — (ieft) Relationship between flare scale height and column density in the lower 
atmosphere of II Peg, using the plasma temperature and equivalent widths appropriate to 
each time slice (values listed in Table 3), for five time intervals in orbits 1 and 2. Above 
column densities of Nh ^LSxlO^"^ cm~^, significant Compton scattering will occur and 6.4 
keV emission will be reduced or eliminated; this region is indicated by the cross-hatching. See 
§4.5.1 for discussion, (right) Predicted equivalent width as a function of N*, the atmospheric 
column density where Ka emission originates, for coUisional ionization formation of the 6.4 
keV Ka line. Top panel shows results for thick-target bremsstrahlung model in Orbit 1; 
bottom panel shows results for Orbit 2. Dashed lines indicate equivalent width measurements 
from Table 2. Crosses indicate values for best-fit 6, Fq from spectral fitting (Table 2), as 
well as maximum and minimum values using uncertainties in S and Fq. 



